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1. INTRODUCTION 

Coupled Map Lattices (CML) have been proposed as simplest models for space-time 
complexity. 

Spatially extended dynamical systems are typically present in many experimental 
situations. 

Examples of such systems can be found in hydrodynamics, plasma physics, chem- 
ical reactions or biological systems. For those non-linear systems, CML gives a useful 
alternative point of view to the more conventional tool of partial differential equations. 

The advantage of CML is based on the use of local maps with well controlled dynam- 
ical properties, that are synchronously updated and with suitably chosen local or global 
coupling [1-4]. 

Different models of CML have been widely investigated from various point of view: 
existence and properties of an invariant measure for infinite lattices at low coupling [5], 
stability and universal properties of homogeneous states [3], [7], [8], [9], study of the 
dynamical regime by means of a bi-orthogonal decomposition [10], [11], [17] and, probably 
the most important issue in view of the physical applications as they appear, for instance in 
Rayleigh-Benard convection [12], the understanding of spatiotemporal intermittency [1-3] 
and [13-16]. The emergence of spatial coherent structures from a local chaotic dynamics 
in CML was discussed in [6]. On the other hand, the coexistence of different domains in 
the lattice with coherent dynamical behaviour, separated by boundaries (kinks), was one 
of the first and fascinating features observed in these models. 

The relation of this situation with spatiotemporal intermittency has been mentioned 
by several authors [1], [4], [15] and was named "natural wavelengths" by the later authors. 

In a way, this mechanism can be described as follows: the size of the different coherent 
regions is modified by a slow displacement of the kinks, until a destabilization of the 
periodic structure occurs inside them. Then, after a transient regime of variable duration, 
a new periodic state is restored in those regions. In fig. 1-a, b and c we present three plots 
of the dynamics for typical situations. 

There, we can see spatiotemporal representations of a one-dimensional lattice - see 
(2.1) for an exact definition of the model - for different values of the parameters. The stro- 
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boscopic period AT = 512 or 256 is chosen according to the time scale of the modifications 
observed in the lattice. 

As already pointed out by Kaneko [1], [4] and Keller and Farmer [15], the slow motion 
of the kinks is due to the presence, between the two regions, of a nodal point Xi, lying in 
the neighbourhood of the fixed point of the logistic map. 

In any case, the spatial wavelength is modulated by the motion of the domain walls 
until the destabilization of this wavelength occurs as a threshold phenomenon. 

It will be shown later that bifurcations of the temporal period of the dynamics may 
also occur. 

It is also obvious that, inside each region, the wave number, an integer, must be 
constant between two consecutive unstable situations. Therefore, changing the size of the 
region, will involve a continuous change of the spatial wavelength A, of the corresponding 
configuration. 

From this point of view it seems natural to study the stability of these wavelengths. 
They act as a bifurcation parameter in this transition between local order and chaos as 
was first observed in [15], and therefore we guess that it may help in the understanding of 
spatiotemporal intermittency. 

The aim of the present work is to study analyticaly and numericaly some features of 
this mechanism. Let us mention the two main ingredients used in our paper. 

First, we propose a two-dimensional mapping depending on two parameters, whose 
dynamics (cycles, bifurcations) may describe the main features of the whole lattice. In a 
sense that will become clearer later, we can say that the whole CML dynamics is enslaved 
by the one of this simple model. 

The first parameter of the reduced map is related to the nonlinearity of the local map 
in the CML (here the logistic map) describing the dynamics at each site of the lattice. 
The second parameter contains the diffusive constant governing the interaction between 
neighbouring sites together with the space wavelength and, therefore, is used as a control 
parameter. 

The second main ingredient is used in numerical simulations. Instead of trying to 
describe the "real" situation as it is shown in fig.l, we build an "artificial" recipe where 



the whole lattice is made of only one region and the stability of the space wavelength is 
studied by changing the number of sites in the lattice. 

This is done by means of an " adiabatic compression" described in section 3 below and 
such that the reader can easily imagine by having a look at fig. 2. 

The paper is organized as follows: 

In section 2, after a brief presentation of the models of interest, we use the Fourier 
transform of spatial structures to introduce the reduced model that is described in section 
4. Obviously, this step is reminescent of amplitude equations as they are used in weak 
turbulence [18]. 

In section 5 we present the simplest case for which the reduced model is exact. It 
will be come clear that this case corresponds to the situation described in [6]. A way 
of using this reduced model out of the region of stability, by means of a bi-orthogonal 
decomposition was proposed in [17]. 

In section 6 we deal with the stability in presence of a time-period two dynamics and 
in section 7 of period four. 

Section 8 deals with the case of homogeneous dynamics. 

In section 9 we briefiy compare the previous results with numerical computations of 
Lyapunov exponents and section 10 concludes. 



2. THE MODELS AND THEIR FOURIER ANALYSIS. 

In this work we consider a one-dimensional lattice of length L and a nonlinear map / 
of the unit interval into itself. The state of the system at time t in configuration space is 
given by x*, o < x\ <1, 1 < i < L. The new state at time t + 1 is defined in each site 
by the following convex linear combination of the updated values of neighbour sites : 

xl+' ^{l-e)f {xl) + I {/ {x^,) + f {xU) } (2.1) 



that we often will write : 



X 



t+i 



e 



f {xl) + |A/(xD (2.2) 

where A is the discrete Laplacian operator. 

Here e is the coupling between neighbour sites, < e < 1 and / is chosen as the 
logistic map : 

fix) = ^x{l-x), 0<^i<4 (2.3) 

We always consider periodic boundary conditions in (2.1), namely xl+i = a^i, 

XL = Xq. 

In the following we shall see that most of the results obtained for the model (2.1) are 
still valid for a more general one, where the coupling between sites is given by a diffusion 
operator. This case can be considered as more realistic if one has in mind a limiting process 
to obtain a continuous configuration space. 

For this model, already introduced in [15], we replace formula (2.1) by the following : 

N 

^l^' = T. P^f W+.) (2.4) 

j = -N 

where A^ is the range of the interaction and the pj are real positive numbers satisfying : 

^p, = l (2.5) 

j 

and 

Pj = P-j (2-6) 
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Clearly the model described by (2.1) is a particular case of this one. 
In the numerical simulations of (2.4) we have chosen a sequence of linearly decreasing 
weights : 



P,=c(^l-^^J, 0<\j\<N + l (2.7) 

where c is given by the normalization relation (2.5). 

Finally the parameter n of the logistic map will always verify n > Hc, where fic = 
3, 5699 ... is the critical value above which the map / has a chaotic dynamics. 

In some cases we will have to consider values of n such that /Uc < /U < /Ui = 3, 677 . . . 
for which the dynamics of / is attracted inside two disjoint intervals Ai and A2 whose 
extreme points are determined by the four first iterated by / of the point x = 1/2. 

Then, every point of Ai (resp A2) is mapped by / in a point of A2 (resp Ai). 

For n > fxi this two intervals overlaps. 

Let's now concentrate on the model (2.1). 

Since we are interested in the study of space periodic configurations, it is natural to 
introduce their discrete Fourier transforms. 

In absence of convective terms, we can restrict ourselves to even configurations, bring- 
ing to the following representation : 



L/2 
X^ 

k=0 



S = E4co.(^') (2.8) 



from where immediatly follows : 

L/2 



AxS = E24^o^(^ -1 coJH^) (2.9) 



k=0 



Finally, using the explicit form of the logistic map, we find the equations which give 
the dynamics of the Fourier coefficients. 

For the saving of notations we write ak = a^ and ak = a]^ . 



Therefore we get : 



L/2 

ao = f (ao) - - XI ^^^ + "'L/2 

f 1 ''^'"' 

aL/2 = (1 - 2£) /U < aL/2 (1 - 2ao) - ^ J^ «fc «l/2- 
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and for Q = 1, 2, . . . , L/2 - 1 : 

q-l 
2 



il 
ttg (1 - 2ao) - 2 X] '^'t «g-fc- 
A;=l 

^ g i/2-g I 

^ 2^ afc+i/2-g ai/2-fc - 2^ Ofc+g Ofc > (2.10) 



2 

fc=0 k=l 



where 



a{q,e) = 1 ~2 esin^ (^') (2.11) 



We denote by JF this (L/2 + 1) dimensional map. Notice that, for each q, the non hnear 
part of JF is a convolution product which is proportional to the sum of the a^ aj where 
the indices (i, j) belong to the boundary of a rectangle in the L/2 x L/2 square lattice 
of Fourier modes defined by the corners (0, q) , {q, 0) , (L/2, L/2 — q) , (L/2 — g, L/2). 
This rectangle is reduced to one of the diagonals for each of the exceptional terms ao and 
aL/2- 

It is worth to notice that the coupling parameter e enters in JF only by means of 
(2.11). Therefore, the relation (2.11) has the meaning of a rescaling between e and the set 
of wavenumbers q. 

In order to motivate the introduction later on of a reduced model we make now a 
remark about the relation of equations (2.10) with the corresponding dynamics in physical 
space. 

The first equation of (2.10) describes the dynamics of the Xi spatial mean value since 

1 ^ 

ao=<x>= j'^Xi (2.12) 



On the other hand, writing 

Xi^<x> -\-5i (2.13) 

we can rewrite the first equation in (2.10) as well as the dynamics of the variances 5i. 
Thanks to the chosen boundary conditions we may have < A/(x) >= and therefore we 
get 

<x> = f {< x>) - ^ < d^ > 

'S, = fi{l - 2 < X >) (l + ^A^ S, - fi (l + |a) {6f- < 6^ >) (2.14) 

where, again u, u stand for the values of u at time t and t + 1. 

Coming back to (2.10) we also remark that, if at time to w^ start with a configuration 
for which the only non vanishing coefficients are ao and Qq^ {qo < L/4) then at time to + 1 
the only new created harmonics corresponds to a2qo ■ Therefore if qo and e are such that 

a{2qo,e) = (2.15) 

the only non vanishing amplitude for any time remains a^^. 
In other words, in that situation, 

X* = 4 + a*^ cos ^y^ (2.16) 

is an exact solution. This is the case, for instance for 

Ao = — = 6 and e = 2/3 (2.17) 

Qo 

We shall see that such a selection of a spatial wavelength by means of a non linear 
coupling between the spatial mean value and the amplitude of the oscillation is observed 
in much more general situations than the ones strictly obeying (2.15). 

In order to get an insight into this phenomenon it is useful to look for the behaviour 
of a{q, e) as a function of q, for a given e. As an example this function is plotted in fig. 
3 for £ = 0, 667 and L = 100. The behaviour of a near to the point where a = indicates 
we may expect the damping of all a set of harmonics. 



We can consider the case A = 2 (go = L/2) to be very special, since then the system 
evolves according to only the first two equations in (2.10) where aLi2 = aq^ is coupled 
only to the mean value ao- No condition as (2.15) is required for the existence of an exact 
solution given by (2.16). We will come back to this case in section 5. 

To end up this section we briefiy indicate how the previous results should be modified 
if, instead of the simplest model (2.1) we take the more general diffusive model (2.4). 

The dynamics of Fourier coefficients obeys the same equations as (2.10) where we only 
need to replace a (g, e) by : 



TV / ■\ 

(q, M) = l-45^p,sin2(^^j (2.18) 



where A^ is the range of the diffusion operator. Let us notice that (2.18) is a good candidate 

for a limiting process to a continuous space. 

Then, the equations for the dynamics of the mean value and variances, similar to 

(2.14), are the following : 

< x> = /(< x>) - ^ < d'^ > 

(2.19) 
6i = ^{1-2 <x >)V{6^) - ijV{6f- < 6^ >) 

where the diffusion operator V is defined as 

N 
j=-N 

and we used the normalization (2.5) and symmetry (2.6) for the pj as well as the periodic 
boundary conditions. 

Fig. 4 shows a plot of the function a{q, {p}) where the special form of the pj given 
by (2.7) is considered. Compared with Fig. 3, the faster decay of ct as a function of the 
wave number q indicates a selection of larger wavelengths in this case and this is what we 
can observe numerically. 



3. ADIABATIC COMPRESSION : A NUMERICAL RECIPE 

As it was alredy mentioned in the introduction about the "natural" model, (i.e. the 
one build on a lattice of constant size L) , the spatial wavelengths which are present during 
the dynamics are modulated by the motion of the their respective domain walls. Our aim is 
to study the instabilities which give rise to these spatiotemporal oscillations as the source 
of spatiotemporal intermittency. 

In order to compare this situation in the lattice with the behaviour of the simplified 
model presented in Section 4 we propose to isolate that phenomenon by means of a nu- 
merical experiment. We call it adiabatic compression of the lattice that we will describe 
now. 

When the lattice of size L reaches a periodic state (eventually after a transient regime) 
we carry out the following transformation : 

(i) L becomes L' — L — 1 

(ii) Xi becomes x'^ = Xi + {i — l){xi+i — Xi)/{L' — 1) 

(iii) the periodic boundary conditions are restored for the new lattice. (3.1) 

Notice that this transformation preserves the wavenumber but may change the corre- 
sponding wavelength. 

After relaxation we may observe if the new wavelength of spatial oscillations is stable 
or not. 

We repeat the operation until we notice a transition to a disordered state. Then the 
critical value Xcr is given by Xcr — Li/q where q is the wavenumber and Li the actual 
length of the lattice. 

After this transition the lattice will relax (with constant length Li) to a new periodic 
state with new wavenumber qi and, therefore, a wavelength Ai = Li/qi. The same process 
may be repeated now. 

Coming back to fig. 2. a we can see there, two adiabatic compressions followed by the 
correspondent instability transitions. 
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In this example the model is as (2.1) with // = 3.63 and e = 0.667. The initial number 
of sites in the lattice is L = 85 and the wavenumber is g = 10 (A = 8.5). After adiabatic 
compression the first destabilization occurs for Li = 82 and then Acr = 8.2. 

After relaxation the new wavenumber is gi = 7 and therefore A = 82/7 = 11.7. 

After a new adiabatic compression the system undergoes another destabilization for 
L2 = 57 and therefore A'^^ = 57/7 = 8.1 which now relaxes to a new wavenumber ^2 = 5 
which corresponds to A' = 11.4. 

Fig. 2.b gives more details on the first described adiabatic compression, since we use 
now a stroboscopic period AT = 32 instead of the AT = 256 used in fig. 2. a. We can 
see how the initial perturbation propagates along the lattice and then relaxes to the new 
laminar state with qi = 7. 

One may observe that both destabilizations are obtained for close values of 
A (Acr = 8.2 and A^^ = 8.1) as well as are close the two values of the spatial wavelengths 
selected after relaxation (A = 11.7 and A' = 11.4). 

Notice that (Acr — Acr) ^s well as (A — A') are as small as they are allowed to by the 
rational approximation determined by the maximum size of the lattice in this case and the 
correspondent integer wavenumbers. So we may wonder if by taking a much larger lattice 
at the begining a better agreement of these values may be found. This is effectively true, 
but the price to pay is an extremely larger relaxation time (transients) and so even for a 
small gain of precision, we refrain to pursue in this direction. 

Instead we may notice the time period of the dynamics represented in fig. 2 also 
undergoes a bifurcation above A = 11, a fact that we will adress later on. The reason why 
this is not observed in fig. 2 is simply because there, the stroboscopic periods are always 
multiples of the time periods of the dynamics. 

We will use several times the adiabatic compression to test the predictions of the 
simplified model described in the following section, model where A is taken as the varying 
control parameter. 
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4. A REDUCED MODEL 

We define the following two-dimensional dynamical system : 

(4.1) 



Xt+i = f{X,)-^Y^ 



y^+i = a{\e)f'{Xt)Yt 

where a (A, e) = 1 — 2e sin^ (tt/A) and / is defined as the logistic map with parameter 
/U as before. Let's denote F the map defined by (4.1). Remark that now a is considered 
as a function of continuous parameters A and e. 

The definition (4.1) comes from a matching of equations (2.10) describing the dynamics 
of the Fourier modes and the corresponding formulation (2.14) in terms of mean and 
variances. 

It can also be deduced from (2.10) by means of the first terms of an expansion for 
small amplitudes compared with the mean value, in which case Xt represents the space 
mean value of the configuration at time t and Yt the global amplitude of the variance at 
the same time. 

Even if the real meaning of Yt is to be taken with caution, we got a good agreement 
of the values of Yj according the map (4.1) and the variance of the configuration at same 
time as it can be computed form direct numerical simulation with the whole lattice. 

The only exception is for A = 2 where the factor |u/2 in the first equation of (4.1) 
must be replaced by ^ if we want to maintain this interpretation for Y . This is due to the 
periodic boundary conditions. Instead of changing the form of the map F for this special 
case A = 2, we prefer to keep it in the same form as (4.1) since this simply corresponds to 
a rescaling Y' = Y/\^ which obviously do not change the properties of (4.1). 

Any way, we take (4.1) as a starting point and we will show how it behaves on the 
parameter space. Different bifurcations and dynamical regimes of (4.1) will be interpreted 
and then compared to the corresponding ones for the whole lattice. 

Notice the map F depends on three parameters : /U, e and A, but the later two enter 
only in (4.1) by means of the function a{X,e). Therefore F really depends only on two 
parameters : fx and a. 
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We shall compute different cycles of F (for different temporal periods T) and determine 
their stability domains in parameter space. 

When expressed as a function of e the result will be compared to numerical simulations 
showing bifurcation diagrams of the whole CML as a function of e. 

On the other hand, as already written, when the results are expressed as a function 
of A, they will be confirmed by applying to the CML the adiabatic compression defined in 
Section 3. 

Before going in the more detailed study of the various dynamical regimes of the map 
F, we notice a special feature of it that helps to understand some of the results described 
below. 

The second equation in (4.1) only fixes the ratio Y^+i/^t and terefore is invariant 
by scaling of the variable Y as far as the condition < y < 1 is fulfilled. This allows 
the translation described by the first equation to stabilize the logistic map governing the 
dynamics of X. 

Then, the coupling between the two variables is closed because, in order that Yt stays 
in a bounded domain, this stable orbit for X should be such that af'{Xt) oscillates around 
1 (the case where it is always, in modulus, strictly less than 1 is addressed in Section 8). 

Therefore, at the end we see that the dynamics is governed by a balance of the size 
of the oscillations of the mean variable X and the one of the variance amplitude that may 
stabilize each other. 
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5. THE SIMPLEST CASE. 

We begin our analysis with the simplest dynamics of map (4.1), namely fixed points 
(T = l). 

The coordinates of such a fixed point obey the following system of equations : 

X = f{X)-^Y^ 

(5.1) 
Y = af{X).Y 

From the second equation of (5.1) we can see that an obvious fixed point holds with 
Y = 0. But then we discover from the first equation that it corresponds to the (unstable) 
fixed point of the logistic map, X = f{X). According to the previous section, this should 
be an homogeneous state (with vanishing variance) and it is easy to show it is not stable. 
In section 8, we will come back to this issue. 

Non obvious solutions (Y ^ 0) are easily computed as : 



2 \ a[i 



2(U \ Q;/ V <^A* 



(5.2) 

1/2 ^ ' 



Here ct = a (A, e) = 1 — 2£sin^ (^). Notice we are only interested in fixed points for 
which < A < 1 and Y^ < 1. These conditions fix the boundaries, in parameter space, 
for the existence of solutions with T = 1. 

The conditions for the stability of such fixed points are found by computing the 
eigenvalues of the tangent map of (4.1) evaluated on the solutions given by (5.2) : 



Qi 



l/a —liY 

-2a^Y 1 



(5.3) 



As it can be easily shown, for the values of n that we are interested in and unless 
A = 2 or A = 3, the stability of these fixed points will imply e > 1. 
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We have numerically confirmed the predictions also in this case, but we refrain to 
present them here since they go beyond the frame of the present work, see [19]. 

We then focus on the simplest case, A = 2, in order to show we can recover the results 
shown in [6]. 



First we observe that : 



a (2, £) = l-2£ (5.4) 



and by changing variables : 



Y 

B = X j= 

V2 



(5.5) 



(5.6) 



the equations of the fixed points of (4.1) read : 

A = {l-e)f{A)+ef{B) 

B = {l-e)f{B)+ef{A) 

so that (5.1) reduces to a well known form. A and B are the values of the configuration 
for even and odd sites on the lattice and, therefore, the existence of a solution for (5.6) is 
equivalent to the one of the given type for the whole lattice. 

In other words, as already mentioned in section 3, the reduced map (4.1) is exact 
when A = 2. 

From (5.6) we easily find that A and B are the roots of the equation : 

u^ -Su + P = (5.7) 



where 



and 



S = A + B = 1 -^ (5.8) 



P = AB= , ^ , (l -^ ) (5.9) 

^(1-2£)V ^(l-2£).' ^ ^ 



which obviously agree with (5.1) where a is taken as in (5.4). 
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On the other hand, the stabihty of such fixed points is given by the condition 

|z/±|<l (5.10) 

where u± and the two roots of the usual equation deduced from the jacobian of (5.6) 
namely : 

iy^-Tu + D = (5.11) 

where : 

T = {l-e){f{A) + f{B)} 

(5.12) 

D ^ il-2e){f{A).f'{B)} 

with f'{x) = fi{l-2x). 

In fig. 5 we show the diagram of bifurcations in parameter space {p, e), where the 
region corresponding to stable fixed points is filled in grey. 

We also mention it is possible to show that the stability of the solution of (5.6) implies 
the stability of the corresponding solution for the whole lattice (the proof is given in the 
Appendix at the end of the paper) . 

For that purpose, we notice that the eigenvalue equations for the tangent map of the 
CML with a lattice of even length, calculated at the fixed point defined above may be 
reduced, by symmetry, to a set of L/2 eigenvalue problems for 2x2 matrices. Finally, 
inspecting the different cases, we found that, inside the stability domain of the reduced 
system all eigenvalues of the CML tangent map are, in modulus, strictly less than 1. 
Furthermore, the two bifurcations of (5.6) are of different nature. For the lower end point 
of the stability interval the eigenvalues are real and they are complex for the upper end 
point. 

Let us add that, in principle, it is possible to follow the same lines of equations (5.6) 
for solutions of the CML with T = 1 and A > 2. This leads to a set of A-equations replacing 
(5.6) that was studied in [19], but the system became rapidly cumbersome and it can only 
be studied by numerical computation. 

The advantage of our reduced model (4.1) comes from the fact that we always have a 
two-dimensional system, no matter of the wavelength we consider. 



15 



6. TIME-PERIOD TWO DYNAMICS 



We come back again to the reduced model F defined in (4.1) to analyse the case of 
solutions with time period 2. 

Numerical simulation shows the existence of period 2 cycles satisfying the following 
symmetry relation : 



Xt+i 



Yt+i 



Xt 



-Y, 



(6.1) 



No other cycle of period 2 was found for F. 

The coordinates X, Y of such cycles obey then the equations : 

X = /(x)-|y2 

Y = -af{X).Y 

where, again, the dependance of a = 1 — 2£sin^ (^) in A and e has been omitted. 
The solution of (6.2) is immediately found : 



(6.2) 



2 \ an J 
2/j \ ct/ V ctA* 



1/2 



(6.3) 



The corresponding tangent map is, according to (4.1) given by : 



DF{X, Y) 



-fix) -^Y 

[a/"(X).y anX)\ 
and therefore, for Qt = Q2 = DF (X, -Y) . DF {X, Y), one has 



(6.4) 



Q. 



^-2«//2y2 _^h_ljy 



a^ 



a ^ 

2v2 



-2// (1 - a) y 1 - 2a//^y 
and then the determinant and trace of Qt are : 



(6.5) 
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1 ^' 

a 



1 + ^ - Aai?Y^ 



(6.6) 



from which we can compute the eigenvalues z/± as roots of the corresponding usual equation. 

As an example of the result accuracy comparing to the behaviour of the CML, we first 
show a case where ^ = 3.8 and A = 6 are keeped fix and e varies. For this case fig. 6-a 
shows the largest value of the modulus of the eigenvalues of Q2 as a function of e. This 
function is equal to 1 for si = 0.39 .. . in the domain where the eigenvalues are complex 
and for £2 = 0.695 . . . where they are real. 

Fig. 6-b shows the bifurcation diagram of a CML of 60 sites for which the two 
corresponding values of bifurcation are e'l = 0.4 . . . and £2 — 0.684. As explained in the 
introduction larger lattice gives an even better agreement. 

Fig. 6-c gives, for the same interval in e, the bifurcation diagram of the second 
coordinate — y — of the two-dimensional map F, where we can see the relation between 
F and the dynamics of the CML is deeper than the only coincidence of the two bifurcation 
points. This relation will be treated elsewhere using the method described in [17]. 

As a second test we take e = 0.667 to be fixed as well as /U = 3.8, the same value as 
before and taking now A as order parameter. 

The largest value of the modulus of the eigenvalues of Q2 is now plotted as a function 
of A in fig. 7. 

We found Ai = 5.96 and A2 = 7.98 for the corresponding values of bifurcation. Notice 
that, due to the form of the function a (A, e) , Ai corresponds to £2 and A2 to £1. 

Following the discussion of section 3 we compare now these values with the actual 
CML bifurcation diagram by using an adiabatic compression. 

We begin with a snapshot of a lattice which length is L = 63 and a wavenumber q = 9 
(A = 7) as shown in fig. 8-a. Fig. 8-b shows the lattice after adiabatic compression until 
L = 55, where we can see the relaxed state with the same wavenumber is still stable. The 
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destabilization occurs for L = 54 (A'^^ = 6) as it is seen in fig. 8-c. Then tfie system, 
relaxes always with L = 54, to a new stable state as it is shown in fig. 8-d after a transient 
regime of about 4000 iterations. Fig. 8-e shows a superposition of two such states at time 
t and t + 2. Notice that in fig. 8-d and 8-e we are concerned with states of temporal period 
4 and a more complex spatial behaviour which is reminescent of the Eckhaus instability in 
Rayleigh-Benard convection [18]. 

Finally, with reference to fig. 8-f and 8-g, applying again the adiabatic compression 
a state with T = 2 is restored for L = 48 (A2 = 8) showing again a good agreement with 
the prediction of the reduced map F. 

The regions of parameter space where the map F has stable periodic (T = 2) dynamics 
are shown in fig. 9 in the plane (/U, e) for different values of A (A = 2, 4, 6, 8) and in fig. 10 
in the plane (A, e) for n = 3.8. 

We have made a large number of numerical simulations for the CML and found they 
were also in good agreement with the predictions of F as they are shown in these two 
pictures. 

We end up this section by adding that, as already noticed, changing £ in (1 — e) in 
equations (5.4) transforms a fixed point in a cycle of period 2. We notice that this is a 
special case of a more general transformation that takes equation (5.1) of a fixed point of 
F in equation (6.2) of a cycle of period 2, by changing a in —a. 
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7. TIME-PERIOD FOUR DYNAMICS 

Exactly in the same way of section 6 we can study the solutions with period 4 of F 
and then to compare the dynamics of F with CML one. 

Again we observed that the cycles have the symmetry Y -^ —Y of the map F, so we 
only need to describe the cycle four positive numbers Ai, A2, Yi, I2, cycle that is runned 
in the following way : 

(Ai, rO ^ (A2, -Y2) ^ (Ai, -Y,) ^ (A2, Y2) ^ (Ai, rO (7.1) 

This combines the permutation Ai — * A2 between Ai and A2 for the logistic map described 
in section 1 and the flip of period T = 2 studied before. 
Therefore a cycle T = 4 is solution of 

F^{X, Y) = F{X, -Y) (7.2) 

or, in other words, solution of the system : 

^ = / (/w - f ^') - f («• rw- n' 

(7.3) 

«'• /' {fix) - |y2) . /'(A) = -1 

This leads, after straightforward calculations, to the following equation, expressed in 
the variable Z = a. /'(A) : 

Z^ {Z^ + aV (2 - /u) {Z^ - 1) + 4aZ) -1 = (7.4) 

Denoting Zi, Z2 the two real roots of (7.4) we finally get : 



2 V cta , 

1/2 (7-5) 

We observe that, taking Z = ±1 we recover the equivalent expressions for T = 1 and 
T = 2, which is naturally related with the symmetry of the cycles and the definition of Z 
as the ratio of two successive values of Y. 
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The stability of such cycles is studied as before. 

As in the previous section we only give two examples of application since all the other 
cases we have performed show the same accuracy. 

In both cases we take ^ = 3.63, a value for which the attractor of / lies inside two 
disjoint intervals Ai, A2 as mentioned before. 

First we take A = 8 fixed and £ as a parameter. The stability interval of the cycle of 
F, as it is shown in fig. 11-a, has limit points £1 = 0.42 and e = 0.695 whereas the values 
obtained from the bifurcation diagram of the CML, shown in fig. 11-b, gave e'l = 0.41 and 
4 = 0.69. 

Then, as in section 6, we fix e = 0.667 and use A as a control parameter. Prom the 
stability of F we get stability between Ai = 7.85 and A2 = 10.15 as it is shown in fig. 
12, whereas the adiabatic compression shown in fig. 13 give, for the critical values of the 
wavelengths, A; = 65/8 = 8.1 and A^ = 52/5 = 10.4. 
The same comments as in case T = 2, can be done for T = 4. 

We also notice that, when the eigenvalues of Qt cross the unit circle with complex 
values, the whole lattice undergoes a period doubling bifurcation. 

Finally we want to add that we performed the corresponding numerical simulations 
for the values of the parameter for which the reduced map F has cycles of higher order 
and we also find a good agreement with the corresponding dynamics of the CML. 
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8. HOMOGENEOUS DYNAMICS 

A simple observation about the dynamics of the reduced map F is worth to be noticed. 
Since 

sup|/'(A)|=/x (8.1) 

X 

we see that the condition 

|a|.;U<l (8.2) 

imphes, using the second equation of F that Yt ^ when t ^ oo for all initial conditions 
(^0, Yo). 

Therefore, in this case, the asymptotic dynamics is defined as 

Xt+i = f{Xt) 

(8.3) 
Yt =0 

According to the given interpretation of these variables, such a dynamics may cor- 
respond to a dynamics of the CML for which there is synchronized motion (since 1" = 0). 
The corresponding dynamics is obviously given by the one of the logistic map, accord- 
ing to the first equation. 

This regime is already observed in CML, see [20], but suprisingly, the reduced model 
again gives quite accurate predictions in this case. We will address this question elsewhere. 
Before ending this section, let us add the following comment : condition (8.2) is 
fulfilled if 

'^ > e sm^ hr > (8-4) 



2// VA/ 2^ 

which means that, for all £ > there is always a A (e) such as for A > A (e) (8.4) may not 
be true. But this effect can not be seen if the size of the lattice is not large enough and, 
moreover, only for an infinite lattice we should have no such cutoff. At least in the view of 
the reduced model, we see some care need to be taken when dealing with finite size lattices 
instead of the infinite limit case. 
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9. LYAPUNOV EXPONENTS 

Finally, let us present some results obtained by computing maximal Lyapunov expo- 
nents along orbits initialized in the vicinity of the spatiotemporal periodic states which 
have been studied above. 

The aim is to give a more precise meaning for the terms "chaotic" and "laminar" used 
for a phase of a CML, here simply defined by its positive or negative Lyapunov exponents. 

Again we compare the results obtained from the dynamics of the reduced model F 
with the ones coming from JF together with the explicit formulae of his tangent map. 

Fig. 14-a and b shows a particular typical case for F and JF. Here n = 3.8 and A = 6 
whereas e is taken as control parameter. This is the situation described in section 6, where 
T = 2. 

We may observe, again, a precise connection between the two stability intervals. More- 
over, for a smaller interval (0.33.., 0.40..) we see that the Hopf bifurcation of the map F 
corresponds to the period doubling of JF that we have already noticed. 

Many other numerical computations we have performed, for different values of the 
parameters, confirm this agreement between the two dynamics. 
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10. CONCLUDING REMARKS 

In this work we were concerned with the problem of spatioteniporal intermittency as 
it can be observed in CML. Directed by the idea that one of the mechanisms undersetting 
this phenomenon is a competition of local configurations with different wavelengths we 
have addressed the question of the stability of such wavelengths. 

The essential point of our work came from the possibility of describing many features 
of this bifurcation diagram as well as of the corresponding dynamics of the CML by means 
a two-dimensional map. 

The study of this map allows, in a very simple way, to predict the existence and 
stability of large ordered subsets of the lattice. 

This approach propounds a particular coupling between the spatial mean value and 
the variance of the variables in phase space. It also suggests a particular form of a con- 
trol parameter that combines the diffusion constant and the spatial wavelength of the 
configuration. 

Even if this reduced two-dimensional model gives very accurate predictions for the 
CML dynamics we were not able, in the general case, to achieve a complete proof of the 
correspondance, as it is the case for A = 2. 

Further progress will lead to the problem of the study of the full system T as it 
was described in section 2 and would give a more serious basis to the knowledge of the 
conditions of applicability of the two-dimensional map F . 
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APPENDIX : 

We show there the stabihty of state (5.5) of the whole lattice in the stability domain 
of the variable e of the corresponding state of the reduced map (5.1). 

Due to the fact that the tangent map DJ^ of the whole lattice evaluated on the state 
(5.5) commutes with the operator of cyclic permutations of order two, whose eigenvalues 
are the roots of unity : W/t = e^^, k = 1, 2, ... , L/2 ; the eigenvalues problem for DJ^ 
may be easily reduced to the one of the set of L/2 matrices M^ defined by : 

-{l-e)f'{A) ^/'(^)(l+>^) 

_ef'{A)il + Wk) {l-e)f'{B) 
We remark that M^/2 is equivalent to the Jacobian of the reduced map whose trace 
T and determinant D are given by (5.12), and the traces Tk and determinants D^ of M^ 
may be expressed in terms of T and D : 



Mk 



(Al) 



l-2£ 

where : 



Dk=Pk-D (A,2) 



/3. = 1 + 3^ ■ s^n' (^) (A, 3) 

As shown in fig. 5, for ^ > 3.45, the stability domain of the reduced map extends 
from e > 0.8. For such values the trace T is always negative and (3k satisfye for any k in 
1, . . . , L/2, the relations : 

m<i (A 4) 



-/3fc < 7 < 1 (A 5) 



where : 7 = -/3l/4 = - ^{-ll ■ 

Now, with the following notations 



\u\ = maxdu^l, \u |), 

I z/* I = max{max | ^^^ | , max | z^/T | ) , 
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where z/^j. (resp. z/^) are the eigenvalues of matrices M^ (resp. Mj^/2)- 
and : r = - T, A = r^ - 4D, A^ = r^ - 4PkD 

We are going to consider successively the cases : 
1./ D < and 2./ D > with subcases 

2.1/ A > and 2.2/ A < 0, and show in each of them that, under the hypothesis \u\ < 1, 
one finds : |z/u,| < 1. 



1./ If £) < then : |z/| = ^+V^^+4|-P| ^nd, from : Afc = A + 4(/3fc - 1)\D\ and (A,4) it 

follows : Afc < A for any /c = 1, 2, ... , L/2. 

Let us denote by Ki (resp. K2) the subset of the indices k where Afc > (resp. 

Afc < 0) and consider separately Ki and K2 : 

In the first case : Iz/f^^ I = max { ^"'"^ *' ) < kl- 

keKi \ ^ J ~ 

In the second, where the eigenvalues are complex, one has : li^^'^l = max {-/^kD) < 

\f\D\< \^\- 

2./ If Z> > let us consider separately the cases of real and complex eigenvalues v : 

2.1/ if A = r2-4£>> 0, by use of Afc = A + 4D{1-Pk) and (A.4) one finds : Afc > 

and, consequently : l^"*] = max I V'^ - p^ — j ^ ^j^gj^^ with (A, 5) and the relation D < ^ 

one obtains for |z/^| : 

|z/*| < r i — 2 ) ^'^d, with explicit form of r and 7 as functions of £ : |z^*| < 



l-e 



1 + ^=^ <1- 



2e-l V V^e^, 

1.1 1 for A < one has \v\ = yD < 1, we again consider separately the subset Ky, 
and ^4 of indices k for which Afc > and Afc < 0. 

2.2.1/ if Afc > 0, using again (A.5) : \vf^\ = max ( I+VIE^KD\ < "+Vf^ 
and, as a function of e : jz^^^l < ^zrj-(l + V^) < 1- 



2.2.2/ for k such that Afc < (A; G K^) one has : |z/f4| ^ ^ax (V/9fc-D) < v^ 

fcG_R'4 

|z/| which ends the proof. 
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FIGURES CAPTIONS 

Fig. 1 : Three dimensional views of spatiotemporal evolution of a CML of 100 sites {x 
axis) obtained by 100 successive plots (t axis). Half of the lattice has been initialized in 
Ai (resp. A2) with e = 0.667. a) /U = 3.6 and the plotting time period is T = 512. b) 
fi = 3.62 and T = 512. c) fi = 3.63 and T = 256. 

Fig. 2 : Three dimensional views of adiabatic compressions of a CML with initial length 
L = 85, wavenumber q = 10 {X = 8.5) and initial state : x^ = 0.45 + 0.15 sin {2Tvqi/L) . e = 
0.667; n = 3.63. a) Among 120 successive plots (period T = 256) two destabilizations 
occur for L = 82 and L = 57. b) The same as the begining of a), the destabilization 
observed for L = 82 is shown in much more details thanks to the lower plotting period 
r = 32. 

Fig. 3 : A continuous plot of the function a (q, e) defined in (2.11) versus q for e — 0.667 
and L = 100. 

Fig. 4 : A continuous plot of the function a (g, {p}) defined in (2.18) versus q for L = 100 
and linearly decreasing weights p of range 5. 

Fig. 5 : Stability domain (grey area) of the fixed point of map F in parameter space 
£ X //. {e : horizontal axis; // : vertical axis). 

Fig. 6 : Stability of period two cycles versus e with A = 6 and p = 3.8. a) Plot of 
the largest value of the modulus of the eigenvalues of Q2 given in (6.5) versus e. b) 
Bifurcation diagram of a CML of length L = 60 (Horizontal axis : e) initialized for each 
£ as : x° = B + D cos (27rz/L) where B and D are taken in the vicinity of period two 
cycle (6.3) of the reduced map F. For each e a transient of t = 200 has been died and 100 
iterations of the CML are drawn with a plotting period T = 2. c) Bifurcation diagram of 
reduced map F. For each e, Y coordinate of F is plotted for 200 iterations after a transient 
of 200 has been died. 
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Fig. 7 : Plot of the largest value of the modulus of eigenvalues of Q2 versus A with // = 3.8 
and e = 0.667. 

Fig. 8 : A sequence of transitions between period two states of a CML submitted to 
adiabatic compressions, with e = 0.667, // = 3.8 and a stroboscopic period T = 2. a) 
Initial state of a CML of length L = 63 and wavenumber g = 9 (A = 7). b) After successive 
compressions L = 55. The lattice is still stable but a modulation of spatial wavelength 
may observed, c) Now L = 54, a superposition of about 20 successive states shows how 
the destabilization propagates, d) With L = 54, after a transient of about 4000 iterations 
(from c)) the CML relaxes to a new stable state, strongly modulated in space, e) A 
superposition of d) at time t and the state obtained at time t + 2, which shows, in fact, 
a period T = 4. f), g) A state with period T = 2 and wavenumber q = 6 is restored for 
L = 48. 

Fig. 9 : Stability domains of F, (T = 2) for various values of A (WL = 2, 4, 6, 8) in 
parameter space e x n {e : vertical axis, // : horizontal axis). 

Fig. 10 : Stability domain of F, (T = 2) for // — 3.8 in parameter space e x X {e : 
horizontal axis, A : vertical axis). 

Fig. 11 : With fx = 3.63, A = 8 and T = 4 ; a) Plot of the largest value of the modulus of 
eigenvalues of Q4 versus e. b) Bifurcation diagram of a CML of length L = 64. For each 
e 50 iterations are plotted after a transient of t = 200 has been died. 

Fig. 12 : The same as Fig. 11-a but expressed as a function of spatial wavelength A with 
e = 0.667 and the same value of //. 

Fig. 13 : The same scenario as for Fig. 8.a)-g) with e = 0.667, n = 3.63 and a stroboscopic 
period T = 4. a) Initial state of a CML of length L = 72 {q = 8).b) The lattice just before 
destabilization, L = 66. c) Superposition of about 20 successive states when L = 65. d) 
With L = 65, the periodic state obtained after relaxation, e) The state of d) at time t plus 
the state at time t + 4. f) A spatiotemporal periodic state is restored when L = 52. (q = 5) 
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Fig. 14 : Maximal Lyapunov exponants computed along orbits close, at initial time, the 
cycle with T = 2, A = 6, n = 3.8, versus e : a) For reduced map F ; for each e the orbit 
length is 300 after a transient of t = 100 has been died, b) For map JF with 24 Fourier 
modes ; for each e the orbit length is 300 after a transient of t = 50 has been died. 
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